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. Disturbed spiral galaxies with high rates of star formation pervaded clusters 

I of galaxies just a few billion years ago, but nearby clusters exclude spirals in 

I . favor of ellipticals. "Galaxy harassment" (frequent high speed galaxy encoun- 

^ , I ters) drives the morphological transformation of galaxies in clusters, provides 

CJ ■ fuel for quasars in subluminous hosts and leaves detectable debris arcs. Simu- 

' lated images of harassed galaxies are strikingly similar to the distorted spirals 

^ . in clusters at z ~ 0.4 observed by the Hubble Space Telescope. 

■ Clusters of galaxies are unique cosmological laboratories. There are several hundred 

J> ' galaxies moving at relative velocities up to several thousand km in regions no larger 

. than the distance between the Milky Way and its nearest neighbor, the Andromeda galaxy 

I (M31). Clusters of galaxies have been observed at redshifts up to 2^. By understanding their 

(<~^ . evolution over cosmic times, we probe the geometry of the Universe and the development 

^—t ' of its largest structures. 

lO . Nearby rich galaxy clusters are dominated by elliptical "E" and lenticular "SO" galaxies^, 

' mostly low luminosity dwarfs. Twenty years ago Butcher and Oemler^'"^ discovered that 

. clusters at 2>0.4 have a substantial population of "blue galaxies" seen only as fuzzy blobs 

Qh' in their ground based images. Recent Hubble Space Telescope (HST) images revealed that 

Q . the "fuzzy blue blobs" are low luminosity, often disturbed, spiral galaxies "Sp"^"*^. The 

i-i ' HST imaging teams stress that the disturbed blue galaxies are ubiquitous, but few have 

c/3 . other galaxies nearby^ and there were multiple bursts of star formation spanning up to 

' 2 Gyr^. The dramatic transformation of clusters (shown in Figure 1) occured during a 

^ . "look-back time" of just 4-5 billion years, only a few cluster orbital times. In contrast, the 

' morphological fraction in the field shows far less evolution^. 

' Figure 1: The dramatic evolution of clusters from z ~ 0.4 to today is shown in these 

images of their inner 300h~^ kpc, where h is the Hubble constant in units of 100 km s~~^ 
Mpc^^ . The left panel is the nearby Coma cluster at z = 0.023 (courtesy A. Fruchter, B. 
Moore and C. Steidel). Nearly every object surrounding the two central dominant elliptical 
galaxies is an E or SO. The right panel is an HST image of CL0939 at z = 0.41^. This 
cluster is dominated by spiral galaxies; many appear disturbed yet no other galaxies are 
nearby. While others appear normal, they have enhanced rates of star formation betrayed 
by strong H^ and Oil emission lines. The bright elliptical population is already in place 
at this redshift. The difference in galaxy populations is most dramatic at luminosities more 
than two magnitudes fainter than L*, the characteristic break in the luminosity function^^ . 
Fainter than this luminosity, 90% of galaxies in clusters at z ^ 0.4 are bulgeless "Sd" disk 
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systems, whereas 90% of galaxies within nearby clusters are dwarf ellipticals "dE" or dE/SO 
The resolution of both images is ~ 0.5/i~^ kpc. 

Given a mechanism for distorting galaxies and promoting star formation that operates 
when a spiral first enters a cluster, hierarchical clustering models will naturally enhance 
the number of "Butcher-Oemler clusters" at 2 ~ 0.4^^. Proposed mechanisms include: 
mergers''"^'''-^, compression of gas in the high pressure cluster environment^^'^^ and tidal 
compression by the cluster Each of these scenarios can produce star-bursts, but none 
address morphological evolution or identify the remnants of these distorted blue galaxies 
in present day clusters. By analyzing their HST images, Oemler et al "^^ conclude that 
merging is implausible as the blue galaxy fraction is large and the merging probability is 
low. They observe disturbed spirals throughout the cluster, whereas both ram pressure 
stripping and global tides will only operate efficiently near the cluster's center. 

What mechanism drives star-bursts and rapid morphological evolution throughout a 
cluster of galaxies? Although direct mergers are extremely rare, every galaxy experiences a 
high speed close encounter with a bright galaxy once per Gyr. Here, "close" means within 
50 kiloparsecs (several optical radii) and "bright galaxy" is one as luminous as L^,, the 
characteristic break in the galaxy luminosity function^°. (The Milky Way has a luminosity 
of approximately L^.) We use high resolution numerical simulations of galaxies in clusters 
to examine the resulting damage of these close encounters. To distinguish this from other 
collisional effects such as galaxy mergers and galaxy cannibalism, we refer to these frequent 
encounters as "galaxy harassment". 

1 The cluster environment and galaxy harassment 

We simulate the evolution of a small bulgeless spiral galaxy orbiting within a dense cluster 
modeled on Coma (Figure 1). Figure 2 shows an edge on view of the model galaxy along 
with an inventory of its stars, gas and dark matter. Some of our simulations use smoothed 
particle hydrodynamics^^, to evolve the gas component of the disk at resolutions of 100 - 
500 pc. 

We simulated galaxies on circular and elliptical orbits in smooth cluster potentials before 
examining the effects of harassment. The disk galaxy shows little evolution over 5 Gyrs 
when placed on a 450 kpc circular orbit in a smooth cluster potential. The disk becomes 
bar unstable after the first pericentric passage^^ when placed on an eccentric orbit with 
apocenter and pericenter at 600 kpc and 300 kpc respectively. Thereafter, each time the 
galaxy passes through pericenter, the halo loses a small fraction of its mass but stars and 
gas remain bound. The evolution is far more dramatic after including the other cluster 
galaxies. 

Figure 2: An edge on view of a model Sd galaxy. Blue, red and yellow particles make 
up the dark matter, gas and stars respectively. The exponential disk has a scale length of 2. 5 
kpc and scale height of 200 pc. The disk is constructed with a Toomre "stability" parameter 
Q = 1.5 and is run in isolation for 2 Gyr (as shown) before being set into orbit in the 
cluster. The gaseous disk is initially on cold circular orbits. The dark halo is a spherical 
isothermal with a core radius of one kpc that is tidally truncated at the pericenter of the 
galaxy's orbit. Within 20 kpc the ratio of dark matter to stars to gas is 20:5:1 and the total 
mass is ~ lO^^M©. The graph displays the contribution of each component to the rotational 
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velocity of the disk. Using the relationship found by Tully & Fisher^'^ , a circular velocity of 
^ 160 km s~^ corresponds to a galaxy with a luminosity about L*/5. We performed similar 
simulations using model galaxies with a peak circular velocity of ^ 110 km s~^ . The galaxies 
are placed in a model of the Coma cluster that has a total mass within its virialised radius of 
7 X lO^^/i~^M0. Galaxies within the cluster are drawn from a Schechter luminosity function 
normalised to a typical cluster mass-to-light ratio of250h. This results in a cluster with 950 
galaxies brighter than the Magellanic clouds, but only 31 brighter than . Internal velocity 
dispersions are assigned using the observed correlation with luminosity^^ . The assigned 
masses assume that the dark halos are tidally truncated at their pericentric distance reduced 
by a time averaged loss of 25% owing to harassment^^ . The fraction of the cluster's density 
attached to galaxies varies from zero at its center to nearly unity at its virial radius of 1.5 
Mpc. At 300 h~^ kpc, roughly 15% of the mass is bound to galaxies. The rest of the 
cluster mass is in a smoothly distributed background represented by a fixed analytic potential. 
Further details can be found in Moore et al . 

When we include harassment, the havoc wreaked is determined by the masses of bright 
galaxies in clusters. Galaxies in the field have massive dark halos, but there has been 
speculation that these were stripped from individual galaxies within clusters^^. All galaxies 
are tidally limited by the potential field of the cluster. Over the 5 Gyr evolution of the 
cluster, bright galaxies retain more than half of the mass within their tidal radius (measured 
at the pericenter of their orbit) , the rest being liberated by fast encounters with other bright 
galaxics^^ . 

We have been conservative in our simulations to insure that our results are robust. At 
a fixed mean orbital radius or "guiding center" , the effects of harassment becomes stronger 
as orbits become more elongated. In our cluster model, the mean ratio of a perturbing 
galaxy's apocenter to its pericenter is greater than 10-to-l. A galaxy in a rich cluster with 
a guiding center of 450 kpc will have a typical pericenter (rperi) slightly greater than 150 
kpc. The masses we assign galaxies are ~ 2.8 x lO^-'^(rperi/15OA;pc)(L/L*)^/^M0. As a 
result, their mass-to- light ratios are M/L = A4h'^{rperi/W0kpc){L/L^)-^^'^. The luminous 
parts of elliptical galaxies are observed to have mass-to-light ratios of 12/iMq/L0^^, so the 
perturbing galaxies have modest amounts of dark matter. However, we follow the evolution 
of individual harassed galaxies that have apo/peri ratios of 2 (i.e. apocenter at 600 kpc, 
pericenter at 300 kpc). As a result, this galaxy avoids extremes of the cluster distribution 
and starts with a large dark halo mass since its pericentric distance is not far from its 
guiding center. Both effects serve to underestimate the effects of harassment. 

2 The changing morphologies of harassed galaxies 

In Figure 3, we follow the evolution of a harassed galaxy through its slightly eccentric orbit. 
At each phase, we compare with images of galaxies in clusters. Our initial conditions look 
like a normal spiral galaxy found in the field (Figure 3a). Typically, the first encounters 
create "disturbed barred spirals"' with sharp and dramatic features drawn out from the 
dynamically cold disk (Figure 3b) . Tails of material can be pulled out and distorted by the 
tidal field of the cluster (Figure 3c) . The gas distribution often forms ring structures that 
tumble within the stellar bar (Figure 3d) . 
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Figure 3: Comparisons of synthetic images from our simulations with observations of 
harassed galaxies. In the simulated images, the stellar distribution is smoothed and filtered 
to model the resolution of the observations. The limiting brightness contours are about 27 
magnitudes per square arcsecond. (a) The left panel shows the initial model galaxy viewed 
face and edge on. The size of each image is 40 kpc across, (b) The upper image shows a 
disturbed spiral galaxy taken from the Butcher- Oemler cluster CL1447. The lower image 
shows our model galaxy 150 million years after suffering a single strong encounter that 
has pulled the two tails of material from the disk. At this time, the perturbing galaxy has 
already moved over 200 kpc away, (c) The upper image is NGC 4438, a spiral galaxy near 
the center of the Virgo cluster (from Sandage and Bedke^'^ ). The close companion in this 
image is probably not responsible for the disturbance^^ . The lower image shows a snapshot 
of our model galaxy after a Gyr in the cluster. The tidal tails of material pulled from the 
galaxy have been subsequently tidally distorted, (d) The upper image is a spiral galaxy with a 
prominent ring in the distant rich cluster CL0939. The lower image shows the ring structure 
often observed in harassed galaxies, (e) Our model galaxy after 3 Gyrs of evolution. The 
stellar distribution shown here should be compared with the initial model shown in Figure 3a 
and with the dwarf elliptical galaxies shown in Figure 3f. (f) CCD images of dwarf elliptical 
galaxies in the Coma cluster. The pixel scale is roughly the same physical size as shown in 
Figures 3(a-e). 

The evolution is driven by several close encounters that would drive the multiple star- 
bursts inferred from HST data^. Another observational puzzle has been the ubiquity of 
disturbed galaxies with no sign of current interaction^. None of the images of our model 
galaxy has another cluster galaxy within 50 kpc. Over the course of 3 Gyr, the closest 
approach of another galaxy is more than 30 kpc away. Since the relative velocity of strong 
encounters is 1,500 km s~^, and the velocity impulse internal to the galaxy is only <50 
km s^^, the perturbing galaxy moves 100 kpc by the time the disk's response is visible. 

After several strong encounters, the loss of angular momentum to their own dark halos 
and the perturbing galaxies, combined with impulsive heating, leads to a prolate figure 
supported equally by random motions and rotation. The gas sinks to the very center of 
the galaxy and the stellar distribution is heated to the extent that it closely resembles a 
dwarf elliptical, although some retain very thick stellar disks and would be classed as dwarf 
lenticulars. At this stage in the evolution encounters cease to create sharp distortions and 
fail to remove any more material from the compact remnant. 

The final stellar systems have a large degree of rotational support, surface density profiles 
and shapes that are in good agreement with observations. Figure 3e shows the stellar 
configuration after 3 Gyrs and can be compared with real images of dwarf elliptical galaxies 
taken from the Coma cluster (Figure 3f ) . Note that the final photometric axes of the model 
galaxy are tilted with respect to the initial plane of rotation of the disk. 

Using our simulations, we can identify the present-day remnants of the disturbed spirals 
seen at z ~ 0.4. Below L*, two distinct classes of elliptical galaxies are observed. Low 
luminosity Es with high central surface brightness are a rare extension to the sequence 
of bright ellipticals; the archetype is M32. The most numerous galaxies in clusters are 
in a second class of dwarf ellipticals, also known as dwarf spheroidals (dE/dSph). Their 
exponential surface brightness profiles resemble those of spirals, as does the correlation 
of their low central surface brightnesses with total luminosity. They are faint, at least 3 
magnitudes below and as many as 14 magnitudes if one extrapolates to the faintest 
known galaxies in the Local Group, Draco and Ursa Minor^^'^°. Harassment transforms 
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spirals into this latter class of galaxies. 

The observed stellar populations of dE galaxies implies recent star formation activity 
that can easily be understood in our model as a result of recent encounters with cluster 
galaxies. Harassed Sd spiral galaxies undergo a remarkable transformation from one mor- 
phological class to another without any merging taking place. Their dynamical states can 
account for all of the dissimilarities between dwarf elliptical and normal elliptical galaxies. 
Harassment provides the link between the dominant populations of galaxies in clusters at 
z ~ 0.4 and the present-day. 

3 Other effects of harassment: tidal debris and quasar fuel 

Stellar and gaseous material torn from the disk during violent encounters creates debris 
tails that lead and follow the galaxy's orbit through the cluster. At pericenter, the tidal 
debris creates a giant arc that could be mistaken for a gravitational arc when viewed face- 
on (Figure 4). However, the transverse size is slightly greater than a lensed galaxy and its 
redshift will match the cluster. The harassment of the debris tails should create tidal shocks 
that promote the formation of dwarf galaxies, as seen in observations and simulations of 
tails in merging galaxies^^'^^. Low surface brightness features punctuated with dwarf galaxy 
formation should be detectable in most HST images of z ~ 0.4 clusters. 

Figure 4: The smoothed surface brightness of the stellar tidal debris after 4 Gyrs of 
evolution. The image is 2 Mpc across and the intensity of the colour shows the logarithm 
of the smoothed stellar surface density plotted wherever the surface brightness mb < 30 
magnitudes per square arcsecond. The white dots are individual particles from the galaxy's 
dark halo. The two long tails were stripped by strong encounters with other cluster members 
coupled with the mean tidal field of the cluster. At various positions along the orbital path 
the stripped stars arrange themselves into long arcs following the galaxy's orbit. HST images 
should show the brightest features (where the colour is brighter than orange). 

Future observations of intra-cluster light coupled with simulations covering the large 
parameter space of orbits and luminosities will provide interesting constraints on our evo- 
lutionary scenario. Current observations are conflicting, anywhere from 0% to 30% of the 
total cluster luminosity could reside in a diffuse component^^'^^. In our cluster model, 20% 
of the light starts off in galaxies fainter than L*/5. After 5 Gyrs, ~20% of the stars are 
lost to the intra-clustcr medium from our small Sd galaxies. The bulk of the evolution is 
driven by the few (< 5) strong encounters with galaxies brighter than ~ L^,. As a result, 
the evolution is chaotic: whereas one fragile disk galaxy can avoid strong encounters for a 
few Gyrs, another may be completely destroyed. The total quantity of stripped material 
can vary from a few percent, to rarer cases where the entire galaxy is disrupted after several 
strong encounters. 

Recent HST images of quasars suggest that many quasar hosts are not galaxies as 
luminous as This is surprising as simple energy considerations imply that quasars 

need IO^Mq - IO^Mq to fuel the black hole engine. To be conservative, one would prefer 
models where no more than 10% of a galaxy's gas must be channeled to the center. This 
argues for hosts at least as large as our own Milky Way, which has a few times 10^ Mq of 
gas. From 8 HST images of low redshifts quasars, Bahcall et al found that 5 of the host 
galaxies must be 0.5 to 1.5 magnitudes fainter than L*. Spiral galaxies of this luminosity 
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have gas masses of order 10 Mq, a quasar needs nearly all of it for fuel! 

Galaxy harassment can give the quasar the fuel it needs. While 10% of the gas is tidally 
stripped from the galaxy, the remaining 90% of the gas (>1O^M0) sinks to the inner few 
hundred parsecs — the resolution of our simulation — by rapidly losing angular momentum 
to the perturbing galaxies, dark halo and stellar bar. 

Clearly, many quasars have luminous hosts^^. Also, quasars are known to avoid rich 
clusters"''^, so how could they be harassed? Oddly, three of the five quasars with subluminous 
hosts found by Bahcall et al lie in clusters. A fourth is probably in a cluster and the 
environment of the fifth has not been studied. Harassed galaxies provide ideal hosts for 
quasars at intermediate redshifts known to lie in subluminous galaxies^^ . 

4 Concluding Remarks 

Galaxies are metamorphised by their mutual interactions. "Merging" of spirals in groups 
creates bright ellipticals^^. In a cluster, one of these "cannablises " its neighbors to become 
the giant central ellpitical^*^ . The dwarf ellipticals are created by the harassment of low 
luminosity spirals. Harassment has the potential to change any internal property of a 
galaxy within a cluster including the gas distribution and content, the orbital distribution 
of stars and the overall shape. Our first examination has only touched on some of the most 
dramatic changes, the phenomenology of harassment promises to be even richer than that 
of merging and cannibalism. 
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